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lit-  this  paper,  me  consider^the  problem  of  estimating  a  survival  curve 
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from  randomly  right  censored  data  when  it  is  known  to  have^e^  Increasing 

Failure  Rate  Average  (IFRA),  or  to  be^fb)'  New  Better  than  Used  (NBU),  Let 
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Frt(t)  be  the  product-limit  estimator  (PL-estimator)  of  Kaplan  and  Meier  for 
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the  life  distribution.  Since  F^(t)  never  has  the  IFRA  property  and  may  not  be 
NBU,  we  modify  F^(t)  to  have  the  desired  IFRA  (NBU)  properties. 


The  modified  estimators  are  easy  to  compute  and  , under  mild  conditions, 

t|>e-  \  ;A  * 

are  shown  to  be  asymptotically  t0  F^(t)  on  compact  intervals. 


Thus  the  modified  estimators  share  the  asymptotic  properties  of  the 
PL-estimator  F^(t). 
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1.  INTRODUCTION  AND  SUMMARY 

In  reliability  theory  and  survival  analysis,  it  is  often  desirable  to 
estimate  survival  curves  or  equivalently,  life  distributions.  In  some 
circumstances  the  lifetime  X^  of  the  i^  item  is  not  observed,  rather  we  only 
know  that  it  exceeded  a  time  Y*.  For  example,  in  clinical  trials  patients  may 
move  away  or  die  of  some  other  causes  and  therefore  are  lost  to  the  study. 

Let  (X^,  Y^),  1*1,  . ...  n  be  i.i.d.  random  variables  with  X^  and  Y^ 

independent  for  each  i.  Let  F(t)  -  P(Xj  <  t)  and  G(t)  -  P(Yi  <  t)  denote  the 

distribution  function  of  Xj  and  Yj  respectively.  F(t)  is  referred  to  as  the 
life  distribution  end  G,  the  censoring  distribution.  In  random  censorship 
model  one  observes  (Z±,  6*),  i-1,  ....  n,  where  Z*  -  Min  (X*,  Y*), 

6^  ■  I(X*  <  Y*)  and,  1(A)  is  the  indicator  function  of  a  set  A.  Techniques 
for  estimating  F(t)  using  (Zj,  6^),  i**l,  ...,  n,  have  been  known  for  a  long 
time  only  recently  has  there  been  much  concern  with  estimating  F(t)  when  it  is 

known  to  belong  to  a  certain  nonparametric  class  of  distributions. 

A  variety  of  such  classes  which  arise  naturally  in  practice  are  given  in 
Barlow  and  Proschan  (1981).  Lo  and  Phadia  (1984)  treat  the  classes  of  convex 
distributions  and  increasing  failure  rate  distributions.  We  shall  consider 
two  nonparametric  classes  in  this  paper:  (1)  the  class  of  distributions  with 
increasing  failure  rate  average  (IFRA)  and  (2)  the  class  of  distributions  with 
the  "New  Better  than  Used  (NBU)"  property. 

For  a  life  distribution  F(t),  let  F(t)  *  l-F(t)  and  H(t)  *  -log  F(t) 
denote  its  survival  and  hazard  function  respectively.  A  distribution  function 
F(t)  with  F(0)  -  0  is  said  to  be  IFRA  if  H(t)  is  starshaped,  that  is;  H(t)/t 
is  a  nondecreasing  function  of  t.  A  distribution  function  F(t)  is  said  to  be 
NBU  if  F(x+y)  <  F(x)  F(y),  or  equivalently  H(x+Y)  >  H(x)  +  H(Y);  that  is, 

H(t)  is  superadditive.  Note  that  an  IFRA  distribution  is  also  NBU. 
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The  class  of  IFRA  distributions  arises  naturally  in  shock  models  (Esary, 
Marshall  and  Proschan  (1973))  and  is  the  smallest  class  containing  the 
exponential  distributions,  closed  under  the  formation  of  coherent  structure 
and  taking  limits  in  distributions  (Birbaum,  Esary  and  Marshall  (1966)).  The 
NBU  concepts  means  that  the  residual  lifetime  of  a  used  item  tends  to  be 
shorter  than  the  lifetime  of  a  new  item. 

The  product-limit  estimator  (PL-estimator)  Fn(t)  due  to  Kaplan  and  Meier 
(1958),  the  most  cotmnonly  used  nonparametric  estimator  of  F,  will  now  be 
defined.  Let  Z(^)  be  the  it*1  order  statistics  from  the  sample  {Z^,l-l,...,n} 
and  6(i)  be  the  corresponding  indicator  function  associated  with  Zqj.  The 
PL-estimator  is ’defined  by 

(,  6 (i) 

n  (^Sr)  .  «  t  <  Z(n) 

z(i)<t 

0  ,  otherwise. 

Large  sample  properties  of  the  PL-estimator  have  been  studied  extensively 
by  Breslow  and  Crowley  (1974),  Csorgo  and  Horwath  (1983),  Gill  (1983)  and  Lo 
and  Singh  (1984)  among  others.  Moreover,  Wellner  (1982)  showed  that  the 
PL-estimator  posseses  desirable  asymptotic  optimality  properties,  e.g. 
mlnlmax.  In  this  paper  we  modify  the  PL-estimator  (1.1)  so  that  is  has  the 
known  property  (IFRA  or  NBU)  of  F(t)  and  remains  close  to  the  original 
PL-estimator. 

We  shall  restrict  the  estimation  problem  only  on  compact  intervals  [0,T], 
where  T  is  any  point  with  F(T)  <  1  and  G(T)  <  1.  Let  ^(t)  ■  -log  (1-Fn(t)) 
be  the  hazard  function  of  Fn(t).  In  Section  2,  we  construct  estimators  (^(D,,) 
of  the  hazard  function  H(t)  under  the  assumption  that  F  is  IFRA  (NBU)  by 
modifying  H^t).  The  modified  estimator  of  F(t)  is  the  distribution  function 


whose  hazard  function  is  Cn(t)  (Dn(t))  respectively.  The  modified  estimators 
Cn  and  Dn  can  be  expressed  in  close  form  (cf.  (2.1)  and  (2.3))  and  are  easy  to 


compute.  Moreover, 

Sup  |  Cn(t)  -  H(t)  |  <  Sup  j  Hn(t)  -  H(t)  |. 

0<t<T  CKt<T 

Hence  Cn(t)  is  closer  to  H(t)  than  Hn(t)  in  the  sense  of  Kolmogorov  distance. 

The  main  results  of  the  paper  are  Theorems  4.1  and  4.2  where  we  show  that 
under  mild  conditions. 

Sup  n^/2  j  (^(t)  -  Hn(t)  j  and  Sup  n*/2  j  ^(t)  -  Hn(t)  |  tend  to  zero  in 
0<t<T  0<t<T 

probability.  This  implies  then  the  asymptotic  behavior  of  our  modified 
estimators  are  the  same  as  that  of  the  PL-estimator. 

The  proofs  of  Theorems  4.1  and  4.2  utilize  an  1.1. d.  representation  of 
Hn(t)  by  lo  and  Singh  (1984).  Relevant  results  are  summarized  in  Section  3. 

2.  DESCRIPTION  OF  THE  ESTIMATORS 

In  this  section  we  shall  modify  the  PL-estimator  (defined  in  (1.1))  so  it 
has  the  desired  IFRA  or  NBU  property.  We  shall  construct  the  estimators  on 
the  interval  10,T], 

First  consider  the  case  when  F  is  known  to  be  IFRA. 

Let  <^(0  -  Sup  (h(t)  :  h(t)  <  ^(t)  for  0  <  t  <  T,  where  h(t)  is 
starshaped  on  [0, T]},  i.e.  C„(t)  is  the  greatest  starshaped  minorant  (GSM)  of 
Hn(t)  on  [0,T],  It  Is  easy  to  check  that  (^(t)  is  starshaped.  A  closed  form 
of  Cn  can  be  obtained  as  follows: 
n 

Let  tn  *  Z  1(6  ■  1,  Z  <  T)  be  the  number  of  uncensored  observations  in  the 

i-1  1  1 
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interval  [0,T]  and  (Z?,  zP,  . ..,  Z°)  be  the  uncensored  observations  from 

i  z  n 

{Zi,  . Zn}  which  correspond  to  actual  lifetimes,  and  are  in  the  interval 
10,T],  The  PL-estlmator  assigns  positive  mass  on  the  Interval  [0,T]  only  to 
those  points  {Z^,  1-1,  ...,  m).  Let  be  the  smallest  slope  of  all  the  lines 

connecting  the  origin  and  the  point  (Z^,  Hn(Z^-)),  for  k- j ,  j+1,  ...,  m+1, 
where , 

Z^j  -  T,  Hj,  is  taken  to  be  right  continuous  and  hence  H^Z^-)  »  H^Z^Pj). 

That  is,  a  -  min  {H^Zp-l/Z?}  . 

J  j<k<mfl  K 

Using  the  results  of  Wang  (1984a)  in  the  i.l.d.  case,  it  is  easy  to  check  that 
the  modified  estimator  Cp,  which  is  the  GSM  of  H„  Is  a  piecewise  linear 
function  of  the  form: 

(2.1)  Cp(t)  -  <ijt  ,  for  2^  <  t  <  Z®  ,  j-1 . mfl. 

Thus  Cn( t )  has  a  close  form  expression  and  is  easy  to  compute.  The  follow¬ 
ing  lemmas  asserts  that  Cp  is  closer  to  any  starshaped  function  on  [0,T]  than 
Hn,  and  hence  is  closer  to  the  true  hazard  function  H(t)  than  Hp. 

Lemma  2.1.  For  any  starshaped  function  $(t)  on  [0,T], 

Sup  |  Cp(t)  -  4(t)  |  <  Sup  |  Hp(t)  -  4(t)  |  . 

0<t<T  0<t<T 

PROOF:  Let  Sup  |  Hn(t)  -  <J»(t)  |  -  A. 

0<t<T 

If  A  »  «,  there  is  nothing  to  be  proved.  If  A  <  <*>,  4(t)  -A  is  starshaped  on 
[0,T]  and  Hp(t)  >  $(t)  -A.  From  the  definition  of  Cp  we  have  Cp(t)  >  $>(t)-A. 
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Since  Cn  is  a  minorant  of  Hjj,  we  have  C„(t)  <  Hn(t)  <  4>(t)  +  A,  so  the  lemma 
holds.  ■ 

Since  R  is  starshaped  we  have 

Corollary  2.1  Sup  |  Cn(t)  -  H(t)  |  <  Sup  |  H^t)  -  H(t)  |. 

0<t<T  (Xt<T 

Thus  Cj,  is  a  better  estimator  of  H  than  the  hazard  function,  of  the 
PL-estlmator.  A  modified  estimator  of  7(t)  is  then  taken  to  be  the 
distribution  Fn(t)  ■  1  -  exp  {-Cn(t)}.  Note  that  Fn  is  not  necessarily  a 
better  estimate  of  F(t)  than  Fn  although  its  hazard  function  is  a  better 
estimate  of  H(t)  than  Hn.  Theorem  4.1  however  guarantees  its  equivalence  to 
the  PL-estlmator,. 

Next,  we  shall  consider  the  case  when  F(t)  is  assumed  to  be  NBU.  Since 

an  NBU  distribution  has  superadditive  hazard  function.  It  is  therefore 

natural  to  consider 

(2.2)  Dn(t)  -  Inf  (H^s+t)  -  H^s)}  , 

0<t<T 

as  our  estimator  of  the  hazard  function  H(t).  The  estimator  Dn  is  analogous 

to  the  estimator  in  Boyles  and  Samaniego  (1984)  developed  for  the  uncensoring 

case.  Their  results  show  that  D„( t)  is  superadditive.  Since  Hn  is  a  step 

function  with  jumps  at  {Z*?,  Z^,  ...,  Z^} ,  to  compute  one  only  needs  to 

1  2  hi 

take  the  inflnum  in  (2.2)  over  these  points,  i.e., 

(2.3)  Dn(t)  -  Inf  {Hn(Z?  +  t)  -  VZ?)},  where  Z°  -  0, 

0<i<m  1 

and  Dq  is  a  step  function  with  jumps  at  points  of  the  form  Z^  -  Z^  for  some 
r  and  s.  Note  that  D^t)  <  Hn(t)  for  all  t,  and  if  Hn  is  superadditive, 

Dn(t)  ■  Hn(t)  for  all  t.  To  estimate  F(t),  we  again  use  the  distribution 
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F*(t)  whose  hazard  function  is  D„,  i.e.,  F  (t)  ■  1  -  exp  (-Dn(t)}. 
n  n  11 

Unlike  the  I FRA  case  we  cannot  show  that  Dn  is  a  better  estimator  of  H 

than  !!„.  However  we  have 

Lemma  2.2 .  For  any  superaddltive  function  $(t)  on  (0,T], 

Sup  1  Dn(t)  -  $(t)  |  <  2  Sup  |  HjjCt)  -  $(t)  j. 

0<t<T  CKt<T 

PROOF;  Let  Sup  |  H^t)  -  $(t)  \  •  A. 

0<t<T 

It  then  follows  that  Dn(t)  ■  Inf  (Hn(s+t)  -  Hn(s)} 

0<s<T 

>  Inf  (♦(s+t)  -  (j>(s)  -  2A} 

0<s<T 

-  $(t)  “  2A 

where  the  last  step  follows  from  the  superadditlvity  of  $(t).  The  lemma  now 
follows  from  the  fact  that 

Dn(t)  <  Hn(t)  <  *(t)  +  A.  a 

Lemma  2.2  implies  that  Uq  is  a  strongly  uniformly  consistent  estimator  of  H(t) 
on  [0,T]  with  at  least  the  same  rate  of  convergence  as  Hn(t).  Theorems  4,  5 
and  6  of  Boyles  and  Samaniego  (1984)  are  immediate  consequence  of  Lemma  2.2. 

3.  PRELIMINARIES 


fcvtSfLV.V. 


As  mentioned  in  Section  1,  the  techniques  of  our  main  results  (Theorems 
4.1  and  4.2)  are  based  on  the  i.i.d.  representations  of  the  PL-process  derived 
by  Lo  and  Singh  (1984).  In  this  section,  we  shall  give  the  results  that  are 
later  needed  in  Section  4  to  establish  the  main  theorems. 


Lemma  3.1.  (Lo  and  Singh  (1984))  Let  ni>  •••,  hn  be  i.i.d.  random  variables 
with  mean  zero,  variance  a 2  and  jml  <  C  for  some  constant  C  for  all  l<Kn, 

For  any  positive  d  and  z  satisfying  Cz  <  d  and  nza^  <  d^  one  has 

n 

P(  1  Z  n  |  <  3d)  <  2  e"2  . 

i-1 

PROOF;  This  is  Lemma  1  of  Lo  and  Singh  (1984).  The  proof  is  based  upon 
Markov's  inequality  and  Taylor’s  expansions.  ■ 

Let  Q(t)  and  Qj(t)  be  the  distribution  (subdistribution)  function  such 
that  Q(t)  *  F(t)  G(t)  and  Oj(t)  -  P(Zj  <  t  and  6^  *  1).  It  is  easy  to 
check  that  Qi(t)  *  G(t)  d  F(t),  and  hence  dQj(t)  -  G(t)dF(t).  Note  that 

Q(t)  is  the  distribution  of  the  observation  Z.  Let  g(t)  -  [0(s)]  ^  dOi(s) 

be  the  asymptotic  variance  of  Hf,  (Breslow  and  Crowley  (1976)).  For  positive 
real  z  and  t,  and  6  taking  values  0  or  1,  let  £  (z,6,t)  -  [0(z)J  1  I  (z<t  and 
6*1)  -  g(zAt).  The  following  lemma,  although  not  stated  explicitely  in  Lo  and 
Singh  (1984),  can  be  proved  by  tracing  the  arguments  of  their  Theorem  1.  The 
proof  can  be  found  in  Proposition  1  of  Lo,  Mack  and  Wang  (1985). 

Lemma  3.2.  If  F  Is  continuous,  for  any  P  >  0,  there  exists  constant  0  >  0 

(depending  on  P)  such  that 

-1  n 

(3.1)  Vt)  “  H(t)  -  n  1  Z  £  (Zit  6t,  t)  +  Rtl(t),  where 

i-1 

(3.2)  P(  Sup  |  Vt)  I  >  9(lon  n/n)3/*  )  -  0(n"p). 

0<t<T 

We  shall  now  assume  that  F(t)  has  bounded  density  function  f(t)  on  [0,T],  and 


let  M  -  Sup  { { f ( t ) |  :  0<t<T}.  Let 

£i(t)  -  S(Zlt  fit,  t),  and  q  -  Q(T)  >  0. 

Lemma  3.3.  (a)  For  any  x  In  [0,T]  and  e  <  (3/2)Mx, 

n  2  2 

(3.3)  P  {  j  n"  I  ^  (x)  |  >  £  }  <  2  exp  (-  n^-  )  . 

(b)  For  any  x  <  y  both  In  {0,T]  and  e  <  (3/4)M(y-x), 

l  n  n  2  2 

(3.4)  P  {  |n  E  lC1(x)  -  Ct(y>]  |  >  e  }  <  2  exp  {-  )  . 

PROOF:  In  order  to  facilitate  the  apllcation  of  Lemma  3.1,  let  tji  -  Ei(x). 
Then  E  tji  ■  0  and 
o2  -  Var  (m)  " 

-  /J  [Q(t)]~2  d  Q^t) 

«  /*  [Q(t)]“2  G(t)  dF(t),  since  dQ^t)  -  G(t)  dF(t), 

<  F(x)q"2 

<  Mxq~ 2  . 

It  can  be  checked  easily  that 

|51(x)|  <  [O(x))"1  +  [Q(x)f2  <  2q"2  .  Taking  C-  2q"2,  d-  (ne)/3  and 

9  7  2  2 

z  -  (ne  q  )/(9Mx),  then  Cz  <  d,  and  nza  -  d  ,  so  Lemma  3.1  applies  to 
n 

|  E  Tj  |  yielding  (3.3). 

1-1 

To  show  (3.4)  let  ni  •  Ci(x)  -  Ej(y).  Again  Et^  -  0  and 
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Var  (t)^)  »  g(x)  +  g(y)  -  2g(xAy) 

-  JJ  I Q<  t )  ] ~ 2  dQ^t) 

-  J%  tQCt)]"2  G(t)  dF(t) 

<  lF(y)  -  F(x) ]  q'2 

<  M(y-x)  q  2. 

When  C  -  4q“2f  d  ■  (ne)/3  and  z  *  (ne^q2)/ (9M(y-x) ] ,  the  conditions  of  Lemma 
3.1  are  satisfied  and  (3.4)  follows.  ■ 

4.  MAIN  RESULTS 

In  this  section  we  shall  show  the  n^/^-equivalence  Df  the  PL-estimator  Fn 

~  * 

and  the  modified  estimator  Fn  and  Fn.  We  shall  assume  that  F(t)  has  bounded 

density  f(t)  on  [0,T],  and  M  -  Sup  Jf(t)|.  Let  ap  *  Inf  (x:  F(x)  >  0}  be 

0<t<T 


the  left  endpoint  of  the  support  of  F,  and  L  «  T  -  ap. 

The  proofs  of  the  main  results  utilize  notions  of  linear  interpolating 
functions  which  we  now  define. 

Let  {1^}  be  a  sequence  of  integers  tending  to  infinity.  For  each  n, 
partition  the  interval  (ap,  Tj  into  kp  equal  length  subintervals 
Uj.  aj+il*  J”°»  1 »  •••»  V*l»  where  a£  -  aF,  a^-  T. 

For  any  function  g  on  [ap,T]  define  its  linear  interpolating  function 
L„g  as 


Lng(«j)  *  g(a“)  for  l»  •••»  *n»  and  LnS^*)  is  linear  on  [a“,  a^] 

for  each  J-0,  1,  ...,  kn~l. 

Thus  Lng  is  a  piecewise  linear  function.  Note  that  :  ^<t<T}  is 

a  stochastic  process  whose  realizations  are  piecewise  linear  functions. 


n  n 


.  *_s  iv  «y.  \ 'A  vX  yiiA .♦>*££/*. '#*V 


WnfattV^WAmAW. 


Let  ^  be  the  event  that  LnHn  Is  starshaped  on  (0,T)  and  Bn  be  the  event 


that  In^  is  superadditive  on  [0(T]. 


We  proceed  with  a  series  of  Leonas. 


Lena  4.1.  If  there  exists  a  constant  X  >  0  such  that  for  any  0  <  x  <  y  <  T, 


(4.1) 


(H(y)/y]  -  [H(x)/x]  >  X(y-x), 


then  for  any  p  >  0,  and  kn  ■  0  ((n/log  n)  3/8)  t 


2  2  3 

l-P(An)  <  4kn  exp  {  -  —  -3-^. -y-  )  +  0(n  P)  ,  for  n  sufficiently  large. 

36MkJ 

n 


PROOF:  Since  l^Hn  is  piecewise  linear,  it  is  starshaped  on  tO,T)  if  and  only 

n  n  ^n 

if  Hn(a")/a"  is 'increasing  in  j.  This  implies  that  1  -  P(A„)  <  Z  P(E.), 

J  J  J«l 


where  Ej  is  the  event  that  Hn(a^)/a”  is  greater  than  Hn(a^"j)/a.j”j  . 

To  compute  P(Ej)  for  fixed  j>l,  let  x  *  ar  Vi  and  hence  y-x  »  L/kj, 


and  y>x>(L/kn).  Now  consider 


P(Ej)  -  P([Hn(y)/y]  <  [H^xVx]  } 

-  P{[Hn(y)  -  H(y)  ]  -  [Hn(x)  -  H(x)]  <  H(x)  -  H(y)  +  V*)  (nO/x) 

<  P{[Hn(y)  -  H(y)J  -  [Hn(x)  -  H(x)J  <  [Hn(x)  -  H(x)]  (y-x)/x  -  Xy(y-x)  ), 


from  (4.1), 


<  P{[H„(y)  -  H(y) }  -  (Hn(x)  -  H(x)J  <  -  (X/2)y(y-x)> 


+  P(lHn(x)  -  H(x) ]  (y-x)/x  >  (X/2)y(y-x)) 


<  P{ |  ~  Z  l£.(y)  -  C,(x)]|  >  (X/4)y(y-x)>  + 
“  1-1  1  1 


P{|Rn(x)  -  Rn(7)|  >  (X/4)y(y-x)>  +  P{  J  i  Z  5  (x)|  >  (X/4)xy) 

n  i-1  1 


+  p{|F„(x)|  >  (X/4)xy  } 

-  n 

<  p{l  i  2  IZA*)  -  Cj(y)  1  j  >  (X/2)  (y-x)2>  +  P{  (Rq(x)  I  >  a/4)(L/k„)2} 

i-1 

+  P{|My>  I  >  (X/4)(L/k„)2  }  +  P{  j  1  Z  lAx)  I  >  (X/2)  x(y-x)} 

n  i-1  1 

+  P{jRn(x)}  >  (X/2)  (L/kn)2),  since  y  >  a£  -  2(y-x)  -  2(L/kn). 

For  k„  «  0  ((n/log  n)  3/8),  Lemma  3.2  implies  that 

(4.2)  P{  Sup  |  Rfi(t)  |  >  (X/4)  (L/kn)2}  -  0(n"P). 

0<t<T 

Since  yx  -  L/kjj  tends  to  zero  as  n  tends  to  «,  we  have  (X/2)  x  (y ~x)  < 
(3/2)Mx  and  (X/2)  (yx)2  <  (3/4)  M  (y-x)  for  n  sufficiently  large. 

I 

Lemma  3.3  now  implies  that 

n  2  2  2 

(4.3)  P{l  {  r  Vx)l  >  <X/2)  x  <*~x>  >  <  2  «P  {  nJ1  S'6MrX>  } 

<  2  exp  {  -  )  ,  and 

36Mk 

n 

(4.4)  P{|  i  E  (C1(x)  -  5t(y)]|  >  (X/2)  (rx)2}  <  2  exp  {  -  ) 

2.  2  3 

<  2  exp  {  -  BA  -X-L  } 

36Mk 

n 

The  lemma  now  follows  from  (4.2),  (4,3)  and  (4.4),  where  we  replaced 
knO(n”p)  by  0(n"P)  since  P  can  be  any  positive  number.  ■ 
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Next  we  obtain  a  similar  result  for  the  event  We  first  find  a 

necessary  and  sufficient  condition  for  the  polygon  I^R^t)  to  be 
superaddltlve  using  the  results  of  Bruckner  (I960,  1962).  Note  that 

{ 0>  {Sj,  j-0,  l,  ...»  kn)  are  the  vertices  of  the  polygon  LftR^t). 

Lemma  4.2.  l^H^t)  Is  superaddltlve  on  (0,T]  If  and  only  if  It  Is 
superaddltlve  on  (a",  J*0,  1,  1^}* 

PROOF:  We  only  need  to  consider  the  case  when  LnR^t)  is  superaddltlve  on 
»  j-0,  1,  •••,  kj,} « 

Let  fi(t)  -  0,  for  0  <  t  <  ay,  and 

f2<t)  -  LnRn(t  +  for  0  <  t  <  L, 

Then  f2(t)  Is  a  polygon  on  (0,L]  with  equally  spaced  vertices,  and 
Theorem  8  of  Bruckner  (1960)  Implies  that  f2(t)  Is  superaddltlve  on  (0,L]. 
Following  the  definitions  of  Bruckner  (1962,  P.  127),  It  can  be  checked  easily 
that  f}  and  f2  form  a  decomposition  of  LnHn(t).  Since  the  minimal 
superaddltlve  extension  of  fj  is  the  zero  function  on  [0,T],  Theorem  1  of 
Bruckner  (1962)  Implies  that  LnRn(t)  Is  superaddltlve  on  [0,T].  ■ 

Lemma  4.3,  If  there  exists  a  constant  a  >  0  such  that  for  any  0<x<y<T, 

(4,5)  H(xfy)  >  H(x)  +  H(y)  +  axy  ,  then  for  any  P  >  0  and 

k„  -  0  ((n/log  n)  3^8), 

2  2  3 

1  -  P  (Bn)  <  4  k^  exp  {-  -n-^  g  ,  ■  }  +  0(n"P)  for  n  sufficiently  large. 

n  72MkJ 

n 


\  J%/’  J . 


..  Cfei'  -itr.  *w.Y 


w‘*  AV^  *1  tr  «3r  *.*>«• 


PROOF:  Lemma  4.2  implies  that, 

1  -  P(Bn)  *  P{  There  exists  i  and  j  such  that 

VaJ)  +  ^(aj)  >  V»J  +  «J>  > 

^n 

<  S  P(E..)  ,  where  E. ,  is  the  event 
i,J>l  3  3 

that  Hn(aJ)  +  Hn(aJ)  >  H^aJ  +  aj)  . 

To  evaluate  P(Ejj)  for  fixed  i  and  j,  let  a”  -  x,  a"  -  y. 

P(Eij)  -  P{Hn(x)  -  H(x)  +  Hn(y)  -  H(y)  >  Hn(x+y)  -  (H(x)  +  H(y)]} 

<  P{Hn(x)  -  H(x)  +  Vy)  “  H(7>  >  V*+y>  “  H(*fy)  +  axy} 

-  PUHnU)  -  H(x)l  -  lHn(x+y)  -  H(x+Y)]  >  -  [^(y)  -  H(y)]  +  axy) 

I 

<  P{ [Hn(x)  -  H(x) ]  -  [Hn(xfy)  -  H(x+y)J  >  <a/2)xy} 

+  PfH^y)  -  H(y)  >  <a/2)xy} 

<  P(  |  7  I  15.  (x)  -  5.(xfy)]|  >  (a/4)(L/kn)y)  +  P{  |  -Z  ?.(y)| 

”  1-1  1  1  n  l 

>  (aMJd/Vy}  +  3  P{  Sup  |R„(t)|  >  (aMKL/kn)2} 

0<t<T 

<  2  «P  (  -  )  +  2  exp  {  .  sa!2!ai!su)!z  ,  +  0(„-K)  , 

by  Lemma  3.3  since  (a/4)(L/kn)y  <  (3/2)  My  <  (3/4)  My  for  large  n, 
2  2.3 

<  4  exp  {  -  -HS-SL —  }  +  0(n”P)  ,  since  y>2(L/k„). 

72Mk 

n 

2 

Lemma  4.3  thus  follows.  Here  we  use  the  same  fact  that  k  0(n"P)  can  be 

n 

replaced  by  0(n~^)  since  the  choice  of  p  can  be  arbitrary.  a 
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The  following  two  Lemmas  are  propositions  2  and  3  of  Vang  (1982). 
Lemma  4.4.  For  any  distribution  function  F  with  hazard  function  H, 


Sup  ul/2  ! 

|  H(t)  -  LjjH(t)  | 

♦  0  Implies  that 

0<t<T 

Sup  n1/2 

1  «n(t)  -  Wt) 

|  ♦  0  in  probability, 

0<t<T 

Lemma  4.5.  Let  H  be  a  differentiable  hazard  function  satisfying 

(4.6)  |H'(x)  -  H’(y)J  <  0  |x-y| ,  for  any  x,  y  In  [0,T),  and  some  constant 

0>O,  then 

Sup  |H(t)  -  LnH(t) |  <  20  (L/kn)2. 

0<t<T 

I 

We  are  now  ready  to  prove  the  main  theorems. 

Theorem  4.1 

Let  F  be  an  IFRA  distribution  function  satisfying  (4.1)  and  (4.6),  then 
1  /2 

Sup  n  |Cn(t)  -  Hn(t)]  0  in  probability. 

0<t<T 

J222£=  ut  ‘n  •  i-iaiAsr;  i  1/3 • 

Setting  P*3  In  Lemma  4.1  Implies  that  l-P(An)  <  n“2  for  n  sufficiently  large. 
When  An  occurs  LqHo  Is  starshaped  so  Lemma  2.1  Implies  that 

sup  |Cn(t)  -  Hn( t) }  <  sup  |  c„(t)  -  i^yt)  |  +  sup  l  i^yt)  -  yt)  | 

0<t<T  (Xt<T  (Xt<T 

<  2  Sup  |  LnHn(t)  -  Hn(t)  |. 

0<t<T 

1/4  -1 

Since  n  y  -*0,  Lemma  4.5  Implies  that 

Sup  n1/2  }  H(t)  -  LnH(t)  |  -*■  0. 

0<t<T 

The  theorem  now  follows  from  Leema  4.4.  ■ 


Corallary  4,1.  Under  the  assumptions  of  Theorem  4.1, 

Sup  |  F  (t)  -  F  (t)  |  ♦  0  in  probability. 

0<t<T  n  n 

PROOF?  Sup  I  P  (t)  -  Fn(t)  I  -  I  exp  (-Cn(t))  -  exp  (-^(t))  1 
0<t<T 

<  |  yt)  -  Vt)  |.  ■ 

Theorem  4.2.  Let  P  be  an  NBU  distribution  function  satisfying  (4.5)  and 
(4.6),  then 

1/2 

Sup  n  |  Dn(t)  -  Hn(t)  J  ♦  0  in  probability. 

0<  t<T 

PROOF:  The  proof  follows  the  derivation  of  Theorem  4.1  utilizing  Lemmas  2.2, 

2  2.3  ... 

4.3,  4.4  and  4.5,  and  choosing  k„  -  h?16M  fog'n^  1  *  • 


Corallary  4.2.  Uhder  the  assumption  of  Theorem  4.2, 

1/2  ,  * 

Sup  n  '  |  F  (t)  -  F  (t)  |  -*■  0  in  probability. 

0<t<T  n  n 

PROOF:  The  proof  is  similar  to  Corallary  4.1. 


Remarks : 

1.  Condition  (4.1)  essentially  means  that  H  is  uniformly  strictly 
starshaped,  or  the  life  distribution  F  is  uniformly  strictly  IFRA.  Condition 
(4.5)  is  a  uniformly  strict  superadditivity  condition  on  the  hazard  function 
H,  which  essentially  means  that  the  life  distribution  F  should  be  uniformlv 
strictly  NBU. 

2.  From  the  proof  of  Theorems  4.1  and  4.2,  it  can  be  seen  that 
1/2 

Sup  n  |  LnHn(t)  -  H^t)  |  also  tends  to  zero  in  probability  for 
0<t<T 


SMW! 


properly  chosen  { k^)  (e.g.  the  choices  of  k„  In  Theorems  4.1.  and  4.2.).  A 
practical  consequence  Is  that.  Instead  of  computing  (^(Df,)  exactly,  one  may  be 
able  to  use  l^H^  as  our  estimator  of  the  hazard  function.  The  estlmtor 
Is  much  easier  to  compute  than  Cn(Dn)  and  has  a  high  probability  (PCA,,)  or 
P(Bn))  of  being  starshaped  (NBU). 

3.  The  techniques  of  this  paper  can  also  be  applied  to  obtain  analogues 
of  Theorems  4.1  and  4.2  for  distributions  with  Decreasing  Failure  Rate  Average 
(DFRA)  or  distributions  which  are  New  Worse  than  Used  (NWU).  Similar  results 
to  Theorems  4.1  and  4.2  can  be  obtained  for  DFRA  and  NBU  distributions 
respectively. 
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